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Point bosons in a one-dimensional box: 
the ground state, excitations and thermodynamics 

Maksim TomchenkcQ 

Bogolyubov Institute for Theoretical Physics, If-b Metrolohichna Str., Kyiv 03680, Ukraine 

(Dated: August 11, 2015) 

We determine the ground-state energy and the effective dispersion law for a one-dimensional 
system of point bosons under zero boundary conditions. The ground-state energy is close to the 
value for a periodic system. But the dispersion law is essentially different from that for a periodic 
system, if the coupling is weak (weak interaction or high concentration) or intermediate. We propose 
also a new method for construction of the thermodynamics for a gas of point bosons. It turns out 
that the difference in the dispersion laws of systems with periodic and zero boundary conditions 
does not lead to a difference in the thermodynamic quantities. In addition, under zero boundary 
conditions, the microscopic sound velocity does not coincide with the macroscopic one. This means 
that either the method of determination of k in the dispersion law E{k) is unsuitable or the low- 
energy excitations are not phonons. 

PACS numbers: 67.25.dt, 67.25.dr 


I. INTRODUCTION II. INITIAL EQUATIONS 


Systems of many particles with contact interaction 
were studied in a lot of works, starting with the work 
by Bethe [I| (see other references in monographs M 
and recent reviews Bi). Models for spinless bosons 
with point interaction were constructed in several main 
works: these are the Girardeau model 0 for impenetra¬ 
ble particles under periodic boundary conditions (bCs), 
the Lieb-Liniger model Q the Lieb model @ for 
penetrable particles under periodic BCs, and Gaudin so¬ 
lutions [13 for zero BCs. Moreover, some results for the 
ground state under zero BCs were obtained by Batchelor 
et al. m- The influence of the boundaries on the ground 
state of a system of point fermions was studied as well 

[Hill. 

In the present work, we investigate the ground-state, 
the dispersion law, and the thermodynamics for spin¬ 
less point bosons under zero BCs. The two following 
results are basic. We give a new method for construc¬ 
tion of the thermodynamics for point bosons and find 
out that the dispersion law under zero BCs differs from 
that under periodic BCs. However, under zero BCs, the 
low-energy elementary excitations are not phonons, if we 
determine k as in the present work. Therefore, the curve 
E{k) under zero BCs is the “dispersion law” in the ef¬ 
fective sense. The nonphononicity arouses the questions, 
which are considered in Section VII. The dispersion re¬ 
lations for SU (2) point bosons under periodic and zero 
BCs were found in [alii- But the results for periodic 
and zero BCs are given at different system parameters; 
therefore, it is hard to see whether dispersion relations 
depend on boundaries at large N, L. 


*Electronic address: mtomchenko@bitp.kiev.ua 


We now recall briefly the basic equations. They are 
presented in the elegant and exact form in the works by 

Gaudin [i[i3- 

Let us consider N bosons, which occupy a one¬ 
dimensional (ID) interval of length L and interact by 
means of a binary repulsive potential in the form of a 
delta-function. The interval can be closed or nonclosed. 
The Schrodinger equation for such system is usually writ¬ 
ten in the form Q 

i ^ i<j 

Under periodic BCs, the solution of this equation for the 
domain xi < X 2 < •.. < xn is the Bethe ansatz 1,113 

N 

'S{k}ixi,...,XN) =^a{P)e 1=^ , ( 2 ) 

p 

where kp^ is equal to one of fci, ..., /cat, and P means 
all permutations of kp The coefficients a{P) were deter¬ 
mined in 0. Under zero BCs, the solution is a superpo¬ 
sition of plane waves @ and all possible reflected waves 

M- 

The energy of a system of point bosons 

E = kI k2 (3) 

Since the energy is completely determined by the values 
of ki, it is sufficient to find the corresponding ki in order 
to obtain the ground-state energy Eq and the dispersion 
law. 

For a per iodic system, the equations for ki take the 
form [T^lldjl 

N h-h 

Lki = 2nIf — 2 ^ arctan — -(4) 

i=i 




2 


If = + i = ( 5 ) 

where Ui are integers. For the ground state, Ui = 0 for 
all i. For < I 2 < ■ ■ ■ < I^, system ([3]) has the unique 
real solution {ki} [l^. Using the equality 

arctana = {tt/ 2)sgn{a) — arctan (1/a), (6) 

Eqs. (H]) can be rewritten in the equivalent form [l^ 

N 

Lki = 2^711 + 2 E ( 7 ) 

k,i Kj 

j=i j 

The solutions of (|3]) and © are the collections of {h}, 
for which kj ^ ki for any j ^ i. Below, we assume 
the ordering ki < k 2 <■■■< k]\[, at which the equality 
kj = —kfuj-i-j holds for the ground state. 

For a swtem with zero BCs, the equations for ki take 
the form |^, 


1. The symmetry-based reasoning indicates that the or¬ 
dering ki < k 2 < ■. ■ < kjq assumes ni < n 2 <■■■< tin. 
Therefore, rii > 1 for all i. The smallest nt are = 1. 
In addition, in the limit c —>■ 0 we should obtain the 
momenta of free particles ki = tt/L, which also requires 
Ui = 1 [l3|. We have studied Eq. ([5]) numerically by 
the Newton method and have found no states with the 
energy, which is lower then the energy of the state with 
Ui = 1 for any i. We will assume that, for any c, the 
ground state corresponds to = 1 for all i. 

1) Ultraweak coupling (c —^ 0). In this case, \c/{ki — 
kj)\ <C 1 for all i ^ j. Therefore, relation (|S]) for the 
ground state changes into 



Making use of the change ki = \fcjL ■ {qi + -k j^fcL)^ we 
obtain the equation 


AT . 

L\ki\ = irn^ + E arctan — 
. 7=1 ^ ' * 


arctan 


\ki\ + \kj 


- \kj\ 

^ \jiii, i = (8) 


where rii are integers, and rii > 1. With regard for ®, 
these equations can be written in the form 


N 


1 


1 


^ - qj Qi + Qj + 27 r/ 

As c —>■ 0, we obtain 

N 


\3^^■ (12) 


1 I (A^— V)yJcL . . /io\ 

^ - Qj 


27r 


L\h\ = ^/f-f^^arctaniEEE 

-I- arctan \ ^ i = 1,..., TV, (9) 

li =ni+i — 1 . ( 10 ) 

The quantities ki are commonly ordered in the following 
way: 0 < |fci| < \k 2 \ < ... < |/cAr|. For brevity, we will 
write ki instead of |fc;|. In what follows, ki means \ki\ 
everywhere under zero BCs. We denote n = N/L and 
^=€171. 


We set qi = qf^'^ + (iV — l)\/^/(27r), then 


(AT) 

q\ = 


N 

E ~(W) 

j=\ ql 




(AT) 


It#*- 


(14) 


As was noticed by Gaudin [T^ , it is the equation for roots 
of the Hermite polynomial H]\[(q). This can be verified 
by substituting the polynomial 


N 

H^{q) = 2^l[{q-qf^) (15) 

1=1 

into the equation for the Hermite polynomials 


III. GROUND STATE OF BOSONS IN A BOX 

For point bosons in a box, the ground state is described 
by Eqs. (151) - (fTUl) with Ui = 1. We can verify that rii = 1 
in the following way (the strict proof is absent, as far as 
we know). Since fci > 0 is the smallest quasimomentum, 
the equation 


HN{q) - 2qHN{q) + 2NHN{q) = 0 (16) 

and setting q = g-. The roots g-of Hermite polyno¬ 
mials satisfy the relations 

N 

(17) 

1=1 


Lki = 7rni -f > I arctan --;-h arctan --— ) 

^ V ki- kj ki + kj J 

is satisfied only for ?ri > 1. Otherwise, the right-hand 
side is negative, but the left-hand side must be posi¬ 
tive. Therefore, the minimum value of ni is equal to 


^' 2 

E(9E) |7V»i=iV(iV-l)/2 + 0(iV-i) (18) 

1=1 

(to obtain the last property, we found numerically 
for N < 100; values of q^^'^ for A < 20 are given in [13 )■ 
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Using these formulae, we find the ground-state energy for 
7V» 1: 


N 


Eo=Y.k^ 

2=1 


A^tt^ 3(A^ — l)cn ^ 
2 


(19) 


where ^ is small. Formula (HU) was previously obtained 
in [l3|. This formula can be obtained also by algebraic 
transformations. In this way, we get 


N 

TTC 

4 r 9 . 


2ki 


L2 


- i ) l'’'- 


( 20 ) 


The numerical solution of Eqs. by the Newton 

method shows that the correction ^ in (HU) can be ne¬ 
glected if 7 <C l/N. For 7 > 1/iV, the quantity Eq is 
close to the Bogolyubov solution for a periodic system: 
Eq fv Ncn. 

2) Weak and intermediate couplings, N~^ ^ 7 ^ 10. 
In [ 1 ^, it was asserted that Eq for a (A^, L, c)-system 
with boundaries coincides with a half of the energy of 
the periodic (2A^, 2L, c)-system with = —1 for i = 
1,... ,N and rii = 1 for i = -|- I,..., 2N. This is 

not quite so, because the first system is described by 
Eq. (fS]), whereas, for the positive quasimomenta of the 
second system, relation © yields, after the appropriate 
reenumeration, the equation 


Lki = > I arctan--;—harctan--— ) Lvi 

^ V ki — kj ki + kj j 

-I- TT-I-arctan, i = l,...,A^, (21) 

which differs from (|8l) by the term arctan (c/2fci). 
Therefore, there is no exact correspondence between 
the (N^, L, c)-system with boundaries and the periodic 
(2N, 2L, c)-system. 

The value of Eq was already calculated in iEini. 
But M. Gaudin and M. Batchelor et al. used a method 
assuming the proximity of the values of Eq under periodic 
and zero BCs. However, these Eq may strongly differ 
from each other. Therefore, we use a method similar to 
that in j^, which does not require the proximity of the 
solution for Eq to that under periodic BCs. 

Under zero BCs, the quasimomenta ki satisfy Eq. ([8]). 
For A^ ^ 1, fci vary smoothly as i increases. Therefore, 
it is convenient to pass in ([5]) and ([5]) from summation to 
integration. Let us consider ki as a function of c = i/N\ 
ki = v{i/N) = u(c). Then 

ki+i -ki = N~^{dv/d<;)\^=i/N = ^/p{h). (22) 

Since dk = dv = (du(c)/dc)dc = {N/p{k))d<;, we have 
dc = p{k)dk/N. This yields 


feiV 1 

J p{k)dk = N j d^ = N, 


ki 


0 


( 23 ) 


1 fciv 

'^2 ^ ^ (24) 

0 fci 

where f{k) is any function, and A<j = 1/A^. This gives 
the normalization for the density of states p{k) and the 
simple rule for the transition from the summation over ki 
to the integration. If we set f{k) = I in ([24|) . we obtain 


With regard to (1^^ . Eq. © can be written in the 
form 


few 

Lki = J p{k)dk ^arctan 


ki-k ^ ki + k\ 

arctan- | 


(25) 


where i = I,..., A^, or 

2 /ci 

Lki ~ arctan-= m 

c 

kN 

ki — k ki + k\ 

— / p{k)dk I arctan-h arctan-| . (26) 

fei 

Considering \ki+i — ki \ small, we subtract the ith equation 
in (HH) from the {i -I- l)th one. In view of (1^^ . we obtain 
A^ — 1 equations for p{ki). They can be written as the 
integral equation 

kN 


^P{d) - j p{k)dk 


+ {k — qY -\- {k + qY 


= L- 


2c 


where q € [fci, /cat] . Since fci is unknown and can be sepa¬ 
rated from 0 by a gap, we write additionally the equation 
for fci, 

few 

fci — fc fci -I- fc' 


Lki ~ ^ ~ J P{k)dk ^arctan ■ 


arctan ■ 


(28) 


where fc 2 = fci -I- I//o(fci). Equations (1^ . (I?7l) . and 
set the complete system of equations for p(k), fci, and 
kN- 

Let us compare them with the equations for the ground 
state of a periodic system. Let N be even. Eor the pos¬ 
itive quasimomenta fci < fc 2 < ... < fcAr /2 of a periodic 
system, the equations read 


'2W/2 


Pp(k)dk = N/2, 


(29) 


fei 


'^N/2 




ki — k 


Lki = ~ J 2pp{k)dk ^arctan ^-h arctan 


ki-\-k 


-f 7 r( 2 i-l), 1 = 1,..., N/2. 


(30) 
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Equation (l30)) follows from (|4]) and ([5]) and yields 


Kn/2 


2Trpp{q) - / 2pp{k)dk [ 


+ 


+ (A: + qy 


\c^ + {k — qY 

= L. 


(31) 


^N/2 


/ k\ — k ki + k 


Lki =— / 2pp{k)dk arctan-h arctan 


fe 2 




„ 2 fci 

+ TT — 2 arctan-, 

c 


(32) 


k 2 = ki + 1/Ppijki). Equations (1^ . (IHTl) . and (15^ form 
the complete system of equations for a periodic system, 
which is written in the form of the equations for a system 
with boundaries. We now make changes 2pp{k) = p(k), 
kN /2 = kN- Then the equations take the form 


kN 


p(k)dk = N, 


(33) 


kN 


TTp{q) — dk 


p{k)c 


p{k)c 


+ {k — qY + {k + qY 


= T(34) 


kN 


I ~/T\ I ki — k ki + k 

Lki = — p{k)dk arctan-h arctan- 

k2 


+ TT — 2 arctan 


2ki 


(35) 


p{k)dk l^arctan 


ki — k 


k 2 = ki + 2/p{ki). Equation (l3^ can be written as 

kN 

Lki = TT + Ai — J 

fei+i/p(fei) 

ki+k\ 

+ arctan- , 


(36) 


. /2fci 1.5 \ 

Ai Ki arctan- 1 -=— 

V c cp(fci)y 

„ 2ki 1.5 

—2 arctan-arctan-=— 

c cp{ki] 


(37) 


terms — and Ai, which are absent, respectively, 

in ([M|) and (E51) . The term enters the com¬ 

bination L — ^ 2 .^ 4 g 2 ■ As N increases, the value of L 
increases as well (at a fixed density n). Therefore, the 
quantity gives an arbitrarily small contribution 

in the limit N ^ oo. The quantity Ai enters the com¬ 
bination Lki — TT — Ai. The numerical analysis shows 
that, in the regime 7 1/JV^ and JV > 1000, the re¬ 

lations cp{ki) ^ 1 and fci/c <C 1 hold, from whence 
I All <C TT (fci can be estimated as ki ^ kN/ 2 /N). In the 
regime 7 1/lV, the relation |Ai| <C Lki is valid as well. 

Therefore, the quantity Ai can be neglected. 

This means that, for 7 ^ 1/lV, the distinction between 
systems (l33l) . (I34l) . (I36l) and (l23l) . (l27l) . (l28l) is negligible. 
Therefore, ki = ki, kjq = ^Ar/ 2 , and p{k) = p{k) = 
2pp{k). Thus, for a system with zero BCs, p{k) is 2 
times larger and fciv is the same, as compared with p(k) 
and /cjv for a system with periodic BCs and the same 
c,N, L. The contribution to Eq is given by positive and 
negative ki under periodic BCs, and only by positive ki 
under zero BCs. Therefore, the values of Eq under zero 
and periodic BCs almost coincide. The difference in these 
energies is small (aE ^ Eq/N) and can be determined 
by the method (see also [il|), in which one should 
take the term arctan (c/ 2 A:i) into account. 

We note that the same equations for zero BCs can be 
deduced by starting from ([ 8 ]). 

3) Strong coupling (7 ^ 1). This case corresponds to 
very large c or L. Consider Eq. ®. The limit c —>■ 00 
means large denominators on the right-hand side of (P|). 
whereas L —>■ 00 means small ki, i.e., small numerators 
on the right-hand side of (O. In both cases, the sum 
on the right-hand side can be neglected. As a result, we 
obtain the solution 


ki = iri/L, 


(38) 


N 

Eo = Y^ k^ = Nir^n'^/S. (39) 

i=l 

The same solution for Eq is obtained for the periodic 
system Q. This is the Girardeau limit Q. Formula (IMl) 
with a subsequent correction was obtained previously by 
another method m- 


IV. DISPERSION LAW 


A. Zero boundary conditions 


Equations (1551) . (1551) . (155)) differ from (1551) . (1571) . (1551) only In order to understand the meaning of an elementary 
by two terms. Equations (I57D and Eq. dSH) include the excitation for a system of point bosons with zero BCs, 
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let us write the system of Eqs. m again: 


N . ^ 

Lki = TTUi + y ( arctan --— 

J^i\ 


arctan 


ki “h kj 


i = l,...,iV. (40) 


This is a system of N equations, where rii are inte¬ 
gers. The ground state corresponds to = 1 for all 
i = 1,..., iV. If at least one rii > 1, we have an excited 
state. Equations (HUl) can be compared with the keys of 
a piano. The pressing of the Zth key can be interpreted 
as a generation of an elementary excitation with ni > 1. 
The pressing of the jth key (j ^ 1) means the generation 
of the second elementary excitation. 

By minimally pressing I last keys, we obtain a config¬ 
uration with ni<:N+i-i = 1, ni>N+i-i = 2, i.e., I excita¬ 
tions with smallest rii > 1. In work such a structure 
is associated with a “hole,” a second type of elementary 
excitations. And excitations with rii^N = I^un > 2 are 
called “particle states” @ . We note that the analysis @ 
was executed in another language, by starting from Eqs. 
(g]) written for the difference ki+i — kt. However, the 
properties of excitations are most clearly seen from Eqs. 
(H0|) . The separation of the excitations into holes and 
particles is based on the analogy with a Fermi system 
and can be carried out in the same way for periodic and 
zero BCs. But we consider that, for a Bose system, it is 
more natural to describe all excitations in a unified way. 
Most simply, the excitation can be associated with the 
clicking of a single key. In what follows, we will define 
the elementary excitations namely so. 

We now hnd the dispersion law for an elementary ex¬ 
citation. For the ground state, we have rii = 1 for all 
f, and some ki are the solutions of (|40)) . For an excited 
state, we write Eqs. (HOI) in the form 


Lki = TTUi -f f arctan ^ 

c \ 

+ arctan^- — ] \j^i, i=l,...,N, (41) 

h + kj J 

where Ui^^N = 1, > 1. In this case, 0 < ki < k 2 < 

... < kN- We set 

u}{ki) =LOi = h- ki. (42) 

At the transition to the excited state, only the A^th equa¬ 
tion in (1101) is changed. Therefore, we may expect @ that 
uii<N are small (|a;i<Ar| <C |A:i|) and ujn is not small. The 
solution agrees with this assumption. In this case, the 
quasimomentum and the excitation energy are as follows: 



uj{k)p{k)dk + ojn, 


(43) 


kN-1 


N ^ 

A = ^(y — fc^) Ri -I-2fcArWAr-I-2 / kuj{k)p{k)dk. 

i=l , 

fcl 

(44) 

From Eqs. dm) with the numbers i = 1,..., iV — 1, we 
now subtract corresponding Eqs. (HOI) . In view of the 
smallness of uii<:N and the nonsmallness of lon, we obtain 
the equations 

N-l f . 

T _ j ~ I [ I 

+ {k + kjY + c2 j 

+ f{ki), (45) 


f{ki) = arctan — 
— arctan 


ki — kjy 
c 


+ arctan — 


ki + kM 
c 


, , -arctan-———, (46) 

ki — rCjv ki kj\i 

i = , N — 1. The similar consideration of the Nth 

equations in (|40|) and (IdTI) gives the dependence wjv(nAr), 
which is unnecessary for finding of E{p). 

In the left- and right-hand sides of (1451) . we add the 
term with i = j, transit from summation to integration 
by rule (1241) . and extend the domain of definition of a;(fc) 
to negative k by the rule a;(—fc) = —uj{k). Then relations 
(1451) yield 

kN-1 


uj{q) L - 


2c 


= c 


dk ■ p{k) 


uj{k) — uj{q) 


-4q^J ' J """^(q-A:)2 + c2 

— kN-1 

+ /(?), (47) 

q S [—kN- 1 , kN-i]. We expand (l46l) in the small param¬ 
eter uJi<^N and obtain 

c c 


/(?) 


arctan - 


ujn + kN + q 
c 

arctan-arctan 

kN - q 

w{q)c 


— arctan 


UJN + kN — q 


kN + q 


(48) 


uj{q)c 


(? -|- [ujn + kN — q)^ (? -I- (ujn + kN + q)^ 

The sum of the first four terms in (1481) is larger than two 
last ones by a factor of ~ \uJN/uj{q)\ ^ 1, and the term 
2cl{c^ + 4(7^) in (|47)) is small as compared with L (for 
7 ^ 1/A^)- We neglect these three small terms. 

For a periodic (c, N, L)-system, the relation Q 

kN 

Pp{k) 


2'Kpp{k) — L = 2c J dk ■ 


ki — — kN 


{q — kY + (? 


(49) 


holds, where 2pp{k) = p{k) (see the previous section). 
Using (l22l) . we obtain 


kN- 


— kw-i 


dk 


p{k) 


{q — ky + c^ 


1 


kN 

— kN 

1 


dk 


p{k) 


{q — ky + (? 


c^ + iq-kNy c'^ + {q + kNy' 


( 50 ) 
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Equalities (|4^ and (l50)) allow us to write (ITT)) in the form 


kN-1 

T^g{q)=c j dkj—^^-^ + f{q) + f{q), (51) 

— kN-1 


where g{q) 


/(<?) 


^{q)p{q): and 

cujjq) cujjq) 

c2 + (g-/cjv)2 c2 + (q + fc7v)^’ 


(52) 


f{q) ~ arctan- - -arctan- - - 

(jJn + kn + q LON + kN — q 

c c 

+ arctan--arctan--. (53) 

kN - q kN + q 

In (IMIl . the function f{q) can be neglected, because 
l/(9)l < l/('7)l- Since f{-q) = -f{q), relation IH]) 
yields g{—q) = —g{q)- We note that the same equations 
are obtained, if we determine ki — ki from Eqs. (O and 

m- 

The dispersion law E{p) can be obtained by formulae 
(l4^ and (|44|) with kN-i ~ /cat — l/p{kN), if we consider 
Wat as a free parameter varying from 0 to oo and, for each 
Wat, find the function g{q) from Eqs. (I5T|) and (l53l) . The 
quantities A:at and p{kN) follow from the Lieb-Liniger 
equations for the ground state @ (Eq. (l23ll with p{k) —>• 
Pp{k), ki = —kN, and Eq. Equations (ISTl) and (1551) 
can be easily solved numerically. For the replacement 
of the integral in (1511) by a sum, we must merely make 
sure that the step is sufficiently small, so that at small c 
the inequality jq — /c"| < c will be valid for at least ten 
numerical points A:". 

The solutions for E{p) are presented in Figs. 1-3 in 
comparison with the solutions for a periodic system Q 
and the Bogolyubov law [3 [13 for the point potential 

Eb{p) = y/p'^ + Acnp"^. (54) 
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FIG. 1: Dispersion laws for point bosons for n = 1 and 1) 
7 = c = 0 . 001 : the curves in the cases of zero (open circles) 
and periodic (open triangles) BCs and the Bogolyubov law 
(crosses); 2) 7 = c = 0.1: zero BCs (circles), periodic BCs 
(triangles), and the Bogolyubov law (stars). 


By comparing Fig. 3 with Figs. 1 and 2, we see the 
dependence of E{p) on n for the same 7 . In the limit 
p —>■ 00 , all curves approach the asymptotics E = p^. It 
is seen from the figures that, for 7 = c/n <C 1 , the curve 
for periodic BCs is close to the Bogolyubov law, but the 
curve for zero BCs is noticeably different from it. In 
particular, the effective sound velocity (i?/p)|p_>o for zero 
BCs is larger than that by Bogolyubov by 1.56 times for 
7 = 0.001 and n = 1; 100 and by 1.49 times for 7 = 0.1 
and n = 1 ; 100 (our numerical calculation for 7 = 0 . 01 , 
n = 1 gave for periodic BCs the sound velocity to be 0.97 
of the Bogolyubov one). As 7 increases, the curves for 
periodic and zero BCs approach each other. As 7 » 1, 
they are close to the Girardeau curve Eg (p) = p"^ + 27rnp 
0 . The relation 7^1 means L ^ 1 or c ^ 1. It is easily 
seen from Eqs. m and (nni that, in these cases, h ~ 
7 r/f/L Ri ki for i < iV, which yields u){k) —?> 0. Therefore, 
we obtain from (|551) and (H51) p = ujn, E{p) = p'^ + 2kNP- 
Since kN\-y^oo —>■ 7 rn Q, we have the Girardeau law. So, 


FIG. 2: Dispersion laws for n = 1 and 1) 7 = c = 10: 
the curves in the cases of zero (open circles) and periodic 
(open triangles) BCs and the Bogolyubov law (crosses); 2) 
7 = c = 100: zero BCs (circles), periodic BCs (triangles), 
and the Bogolyubov law (stars). 


this law holds at c = 00 Q and c < 1,7 1 [3 for 

periodic BGs and at 7 1 for zero BCs. 

The solution E{p) for zero BCs can be written in 
the Bogolyubov form dMl) with the replacement c —>■ 
c- 'd(A:, c, n), where i?(A:,c, n) depends weakly on k and 
strongly on c and n. 

We note that, for 7 < 1 and small or intermediate 
E and p, the main contribution to E and p is given by 
small perturbations cOi^N- For large E and p, the main 
contribution to these quantities is given by wat. In other 
words, the excitations are collective for small E and p, 
and are quasi-one-particle for large E and p. For 7 ^ 1 , 
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FIG. 3: Dispersion laws for n = 100. 1) c = 10, 7 = 0.1: 
the solutions in the cases of zero (open circles) and periodic 
(open triangles) BCs and the Bogolyubov law (crosses). 2) 
c = 1000, 7 = 10: zero BCs (circles), periodic BCs (trian¬ 
gles), and the Bogolyubov law (stars). 

FIG. 4: Dependence p(cJjv) (I43II for 7 = 0.001: 1) zero 
BCs and n = 1, c = 0.001 (open circles); 2) zero BCs and 
n = 100, c = 0.1 (crosses); 3) periodic BCs and n = 1, 
c = 0.001 (circles); 4) periodic BCs and n = 100, c = 0.1 
(stars). 


BCs. For a periodic (c, N, L)-system, the equations are 
deduced exactly in the same way as in the case of zero 
BCs. Starting from Eqs. and ([S]) or from the equiva¬ 
lent equation we obtain 


few 

2^gp{q) = 2c J 

ki — — fcjv 


(55) 


Uiq) ^ 2arctan ^^+^^~^ 

— 2arctan^^^^^—(56) 
c 

where gp{q) = oj{q)pp{q). Equations (|i5)) and (Hill re¬ 
main valid if we replace p{q) —^ Pp{<l) and consider 
kj^-i Ki k]s[ — l/pp(kM)- The values of Pp{q) and k^ 
can be obtained from the Lieb-Liniger equations Q for 
the ground state of a periodic system (Eq. (1^ with 
p(k) Pp{k), ki = —kjq, and Eq. ((4^ ). In Figs. 1-3, 
we give the dispersion laws for periodic BCs, which were 
obtained numerically from Eqs. (|i51) . (H^ . (I55|) . and (I5S1) 
with the indicated changes. 

In the derivation of Eqs. dSl]) and HMD, we consid¬ 
ered fci,..., fcjv_i to be the quasimomenta of a system 
of N interacting atoms (like under zero BCs). Eqs. 
(2.18)-(2.20) from Q were obtained within another ap¬ 
proach, where ki,... ,kN-i were considered to be the 
quasimomenta of a system of — 1 interacting atoms. 
In the first (second) approach, at a ;at > 0 we have 
uj{q) > 0 {uj{q) < 0). The advantage of the second 
approach is that the derivation of equations is simpler. 
However, the ground state corresponds to N interacting 
atoms. Therefore, the first approach is slightly more ex¬ 
act. But the results in both approaches are very close. 
The first approach has the advantage that the properties 
E{ujn —>■ 0) —>■ 0,p(a;Ar —>■ 0) —>■ 0 follow directly from 
the input equations (Ii51) . (Hil) . (1551) . and (1551) . whereas 
the analogous properties in the second approach are not 
obvious and require a bulky proof Q. 


the excitations are quasi-one-particle for any E and p, 
even for p,ujn —>■ 0. For all curves, PalJ^-^o —>■ 0 for any 
c and n. In Fig. 4, we show the dependence p{ujn) (H51) 
under zero and periodic BCs. 


B. Periodic boundary conditions 

The dispersion law for periodic BCs was found by 
E.H. Lieb In order to understand the reason for 
the influence of boundaries on the dispersion law, let 
us compare the formulae obtained for zero and periodic 


V. NUMERICAL SOLUTION BY THE 
NEWTON METHOD 

Trying to solve equations of the form ([5]) numerically, 
we found that this can be performed easily (for N < 1000) 
and with a high accuracy within the Newton method. 
This method frequently requires the proximity of a bare 
solution to the exact one. However, for systems o and 
(O, the method converges also with the quick choice of 
a bare solution. The essence of the method is as follows. 
Two nonlinear equations 

fl{ki,k2)=0, /2(fcl,fc2) = 0 


(57) 







the Fermi form [l^ 


can be approximately written in the form 

+ ^^ikf\kf^) • {k2 - k^f^) = 0. (58) 

analogously for / 2 . Setting kj = k^p {j = 

1,2), we get a linear recurrence relation between k^ and 

■ Under certain conditions, the collection {k^} 
converges with increasing I to the exact solution {kj}. 



ig(L) 


FIG. 5: Dependence of Eo/N on the system size L under 
zero (rhombs) and periodic (circles) BCs for n = 1, 7 = 0.01. 
The value of Eq is obtained by direct solving of the systems of 
equations 0 and (1S]| within the Newton method. The dotted 
line shows the Bogolyubov solution Eq/N = cn(l — 4 y^/ 37 r) 
taken as 1 (Eq/N — 1). 


00 

P J dp\n(l + e^'^, (59) 

— 00 


where e{p) is some effective energy that is not the energy 
of a quasiparticle. To study the influence of the bound¬ 
aries, it is necessary to carry on the whole analysis [T^ 
anew for zero BCs. However, the weakly excited state 
of a quantum liquid can be considered as a number of 
excitations and is most simply described in the language 
of excitations [ 2 ^. Moreover, the excitations are usually 
observed rather than atoms. If we separate the exci¬ 
tations into “holes” and “particles”, then the excitations 
do not obey some simple statistics. The “mixed” descrip¬ 
tion was constructed [l^ probably just for this reason. 
Below, we propose a simpler way to construct the ther¬ 
modynamics, which will allow us to introduce excitations 
in a self-consistent manner. 

In what follows, the formulae are valid for zero and 
periodic BCs. The Gibbs canonical distribution implies 
that the free energy of the system reads 211 


F = —ksTlvL 2_^e 


(60) 


where j enumerates all possible states of the system, and 
Ej is the energy of the system in the jth state. Any 
excited state of the system is uniquely determined by 
the set of numbers {rii} = (rq, n 2 ,..., un) in 0 or ([5]). 
For the ground state, we have the set ..., ). 

Let in system © or ([5]) the number of equations with 
ni ^ be much less than N. The analysis by the 
Newton method shows that, in this case, the energy of 
the system is 


By this method, we found Eq and E{k) for periodic 
and zero BCs by means of the direct solving of Eqs. © 
and ©. To within 1%, the solutions coincide with the 
above-obtained ones. We get also the dependence Eq{L) 
for n = const (see Fig. 5). The difference in the values of 
Eq for small L from i?o(L —>■ 00 ) is the finite-size effect. 
The similar effect was previously studied in Em. 

The Newton method is a powerful tool for solving of 
systems of the form ([4]), 0, 0, 0. 

VI. THERMODYNAMICS 


N 

E{{n,})^Eo + Y,£k{h), ( 61 ) 

= E{n^°\... ... ,n(°^) - Eq, (62) 

where li = ni — , and i is the number of the equation 

in system © or 0. Because all i are equivalent, we 
have ei{li) = e{li). The quantity £{li) coincides with the 
energy of a quasiparticle dm). li can take the zero value, 
then e(0) = 0. Then relation (1501) can be written in the 
form 


In the following section, we will see that the excitations 
are not phonons. Therefore, we cannot observe them by 
the scattering of some particles. It is of interest to clar¬ 
ify whether the difference in dispersion curves for zero 
and periodic BCs leads to a difference in the measured 
thermodynamic quantities. C.N. Yang and C.P. Yang 
M constructed the thermodynamic description of point 
bosons, by using the mixed language of atoms and exci¬ 
tations (see also @). In particular, the free energy has 


- E e(li)/kBT ^ ^ 

F = EQ-kBT\n ^ e •=! . (63) 

In sum (l6^ . we should take only different states into 
account [21| . For example, all states of the form Ui^^j = 
,nj = 2 are equivalent and must be considered as a 
one state. It is difficult to determine such a sum dMl) 
if N is finite. However, this can be easily performed 
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in the limit N = c»,L = oo,N/L = n. For infinite 

N, the different states from the set (Zi,...,Zjv) can be 
enumerated by the set of numbers {r]i}, where rji is the 
occupation number for the ^th state, and I has the same 
values as any li. It is also necessary to replace 

N 

'^e{li)\N^oo ( 64 ) 

i=l I 

where the numbers 77 / for each I can take the values rji = 

O , 1 , 2 ,..., 00 . As a result, we have 

F = Ao-ZcsTln^e ‘ , (65) 

im) 

where { 77 ;} is the set 771 , 772 , • (the enumeration 

coincides with that one for I, see below). Let us rewrite 
F in the form 

00 

F = Eo- UbT XI In XI 

i 111=0 

Summing the geometric progression, we obtain finally 

F = Ao + /fcBTXln(l-e-^6)/fcBTy ^g^^ 

i 

This formula describes the free energy of a system of 
noninteracting bosons with zero chemical potential 
and the additional summand Eq. Formula (EZl) indicates 
that the weakly excited state of a system of point bosons 
can be considered as a number of elementary excitations 
satisfying the Bose statistics. 

This corresponds to the symmetry of wave functions. 
Indeed, the permutation of the Zth and jth excitations 
means the permutation of ni and rij in Eq. © or ([ 8 ]). 
This leads only to the permutation of ki in the complete 
collection |fc,:|. In this case, the Lieb-Liniger @ and 
Gaudin [l0| | wave functions are invariable. In addition, 
the system can possess several excitations with identical 
ni- These properties again indicates that the excitations 
are bosons for any 7 and k. By their properties, the 
excitations are similar to Bogolyubov ones [la, [l^ . 

We note that equality ((HTI) holds if the total number 
of excitations <C N, i.e., if T is low. If the number of 
excitations is of the order of magnitude of iV, one needs 
to consider their interaction. 

For zero BCs, the levels are numbered as follows: I can 
take the values 1,2,... corresponding to the quasimo¬ 
menta fcW = (E^i)l(n^<Ar-i=l.n,v=Z-tl)-(EfcOlK-<iv = l)- 

i i 

For periodic BCs, we have I = ±1,±2,..., = 

(E^i)lK<N-l=O.n;vr=0 - (EfcOlK-<^=0), 

i i 

By the Newton method, we determined the first 1000 
levels (£(?), for periodic and zero BCs for N = 
1000,n = 1, and 7 = 0.01; 1. It turns out that ez{2l) = 
ep{l) = £p{-l), £z{21 - 1) = 0.5(£p(/) -I- £p{l - 1)) (the 
indices p and z mean periodic and zero BCs, and we con¬ 
sider £p[l = 0) = 0). These equalities hold with accuracy 


^ 1/N. If we neglect the small difference £p(/) — £p(Z — l), 
we obtain £^(1) = £p(-l), £ 2 ( 2 ) = £p(l), £ 2 ( 3 ) = 
£p(—2), £ 2 ( 4 ) = £p(2), and so on. In other words, the 
energies of levels under periodic and zero BCs coincide. 
Therefore, the free energies ((63l) for the systems under 
periodic and zero BCs are identical. With regard to the 
difference £p(/) — £p(Z — 1), we obtain a surface correction 
^ F/N to the free energy F. 

In (|67l) . we may pass to the integration with respect to 
p or E: 

00 

Fk,Eq+ j dpC(p)ln(l-e "b'^'’), (68) 


/ dEaE)\n(l-e^), (69) 

— CO 

where ()(pi) = 27 r[L(A:^*+^^ — for periodic BCs and 

C(Pi) = 7 r[L(/c 6 +i) _ for zero BCs, pi = hk^^\ 

C{E) = C(p) • dp/dE. For zero BCs, C(“bl) = C(bl)- 
The numerical analysis indicates that, for periodic BCs, 
— A: 6 ) = 27r/L and C(p) = 1- For zero BCs, the 
step fc 6 +i) _ fclO depends on I and differs from the step 
for periodic BCs (which leads to the different dispersion 
law). In this case, CiE) is the same under periodic and 
zero BCs with deviation < 0.3%. Therefore, relation 
(l69l) allows us again to conclude that E is independent 
of the boundaries. 

Formula (l68)) with ({p) = 1 describes the free energy 
of a gas of excitations in He II [1^ (with the additional 
term Eq independent of T). 

It is of interest that the total number of excitations in 
a gas of point particles is at most N, by definition. For a 
gas of nonpoint particles, the excitation is manifested as 
a multiplier of the total wave function [T^, [23l - l^ , and 
the number of multipliers is unbounded. Moreover, the 
system of energy levels of N point particles has no level 
corresponding to IV -I- 1 phonons in a gas of N nonpoint 
particles. That is, some states of a gas of nonpoint par¬ 
ticles have an analog in a gas of point particles, whereas 
another states do not have. Intuitively it seems that a 
real system of 100 Bose particles can hold 150 oscilla¬ 
tory waves (phonons). It is possible that, for finite N, 
the excitations for point particles are defined not quite 
self-consistently. 


VII. NATURE OF EXCITATIONS, OPEN 
QUESTIONS 

We have shown above that, for a weak and intermedi¬ 
ate couplings, the dispersion law for point bosons does de¬ 
pend on the presence of boundaries, whereas the ground- 
state energy is independent of boundaries (to within the 
surface correction ^ Eq/N). Why is it so as regards 
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equations? As was shown in Sec. Ill, the equations for 
the ground state for systems with zero and periodic BCs 
differ from each other only by two terms, which are small 
for large N and L. But the equations for the dispersion 
law differ from each other strongly. Namely, for peri¬ 
odic BCs, the functions fp{q) and gp{q) in Eqs. ([55]) . 
(1551) are constant-sign for all q. But, under zero BCs, 
f{-Q) = -fiq) and g{-q) = -g{q) (see Eqs. ([51]), ([52])). 
Thus, the equations under periodic and zero BCs differ 
by their symmetries. Therefore, the solutions for the dis¬ 
persion laws are also different. 

However, for zero BCs, the low-lying excitations are 
not phonons. In order to see this, let us compare the 
microscopic sound velocity = dE{p)/dp\p^o with 

the macroscopic one = ffdP/dp,P = —dEo/dV 

{p = mn) [^. Under the periodic BCs, they are identi¬ 
cal Bi- Therefore, the excitations with small p can be 
interpreted as phonons. Under zero BCs, the system is 
characterized by the same Eq and but by the dif¬ 
ferent value of being approximately 1.5 times over 

y^acro at 7 <C 1. Hence, under zero BCs, the excitations 
with small p are not phonons. This is true for 7 < 10. 
For 100 < 7 < 00 , the relation « .^macro jjolds, 

and the excitations are almost phonons. For 7 = 00 , we 
have = ymacro^ excitations can be consid¬ 

ered as phonons. Which structure of the wave function 
(WE) under zero BCs should be in order that an excita¬ 
tion be a phonon? The total WFs are not eigenfunctions 
of the operator of total momentum even for 7 = 00 . But 
the WF can contain a multiplier corresponding to two 
counter-propagating waves. In this case, the WF of a 
low-energy state should have the form 

.. .,xn) = (V'(a;i, . ■ . ,XN\k)+'ipixi, ..., a:Ar|-fc))5'o, 

(70) 

where the function -ip^Xi, ... ,X]s[\k) is an eigenfunction 
of the total operator of momentum with the eigenvalue 
k. The structure of (l70l) is phonon-like. If such represen¬ 
tation is possible, then the excited state has the phonon 
structure and is characterized by the quasimomentum k. 
In this case, the relation = ymacro 3 ijould hold. 

But since is the same as for periodic BCs, the dis¬ 

persion law E{k) should coincide with that for a periodic 
system. We do not know whether a representation of the 
form dZOl) exists. 

According to solutions [H, the boundaries of a 
system of nonpoint bosons affect both the dispersion law 
and the ground-state energy. For a ID system of almost 
point bosons with the weak interaction and zero BCs, the 
following dispersion law is found [ 2 ^ : 


E{p) = _|_ 2cnp2. (71) 

For the solutions [1^, the equality = y^acro j^olds; 

this can be verified for the weak coupling. The dispersion 
law dm is characterized by the sound velocity, which is 
•\/2 times less than the Bogolyubov one (see (1541) b How¬ 
ever, for the point bosons in a box, the effective is 


larger than the Bogolyubov Vg. If there exists a contin¬ 
uous transition from a nonpoint interaction to the point 
one, then the noncoincidence of solutions [1^ with those 
for point bosons indicates the incorrectness of either solu¬ 
tions (23 or the solutions of the present work. However, 


we solved the Gaudin equations for by two different 
methods and are sure in the validity of the solutions. In 
addition, if there is no phonon representation for exci¬ 
tations of a gas of point bosons under zero BCs, we are 
faced with the difficulty for the theory of point bosons. 
Indeed, for the real quantum Bose liquids, the low-lying 
excitations are phonons. This is testified, for example, 
by the experiments on the scattering of neutrons and by 
the measurements of the heat capacity of 3D He H in a 
vessel (zero BCs). The distinction of the one- and three- 
dimensional cases should not be important, because sev¬ 
eral microscopic models of He H [^, l23l - l^ work in ID 
and 3D and give for ID and 3D the solutions of the same 
structure. 


As a possible reason for all disagreements, we can in¬ 
dicate the absence of a continuous transition from a non¬ 
point interaction to the point one, i.e., the anomality of 
the d-function. The d-function is a singular generalized 
function. It is commonly accepted that the replacement 
of a real potential by the d-function is admissible. In par¬ 
ticular, it was proved mathematically [ 2 ^ that the energy 
levels of a system of 3D bosons in a very extended trap 
are close to those of the Lieb-Liniger problem. However, 
the wave functions of N nonpoint and N point particles 
have different forms [ 2 ^. Is it the different forms of the 
same functions or the evidence of the difference of the 
functions? It is known only Q that, under periodic BCs, 
the WFs of point bosons with 7 = 00 can be written as 
the zero approximation for the WFs of nonpoint bosons. 
It is necessary to show that, for an Wparticle ID sys¬ 
tem, all energy levels and the WFs for the almost point 
and point interactions coincide. In the Appendix, this 
is proved for the one-particle problem. The same should 
be proved at least for N = 2 as well. It is of interest 
that, in the 2D- and 3D-spaces, the potential i5(x) has 
no influence on the solutions of the Schrddinger equation 
for some tasks [^ . 

The Muations for spin systems and point bosons are 
similar . Therefore, one can expect that the dispersion 
laws of spin waves under zero and periodic BCs should be 
different (by our method of determination of fc), whereas 
the thermodynamic quantities should coincide. From 
whence, we may conclude that the difference of the curves 
E{k) is unobservable. But it was shown in experiments 
on the scattering of neutrons that the low-lying excita¬ 
tions of magnetics with boundaries (zero BCs) are quite 
observable and have quasimomentum. The possible rea¬ 
son is that the spin wave is accompanied by the sound 
wave. An important point is that the contact Hamil¬ 
tonian describes well the real exchange interaction and 
contains no d-function. Therefore, the solutions should 
correspond to the natural properties, and disagreements 
due to the d-function should not arise. Such a deductive 
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method indicates that the problems of solutions for spin 
waves and point bosons under zero BCs are not, appar¬ 
ently, related to the 5-function. 


VIII. CONCLUSION 

We have obtained two main results: 1) It is found that 
the dispersion law E{k) of a system of point bosons de¬ 
pends strongly on boundaries in the regimes of weak and 
intermediate coupling. 2) The thermodynamics of a gas 
of point bosons is constructed by a new method. Our 
analysis shows that the values of thermodynamic quanti¬ 
ties are independent of the boundaries. By our method of 
determination of the quasimomentum k of an excitation, 
it turns out that, under zero BCs, the low-energy exci¬ 
tations are characterized by a linear dispersion law and 
a nonphonon structure of the wave function. It seems 
strange, because the experiment indicates that the low- 
lying excitations of real uniform quantum liquids with 
zero boundary conditions are phonons. It is possible that 
there exists a way of determination of k under zero BCs, 
for which the low-energy excitations are phonons. Oth¬ 
erwise, the solutions for point bosons do not describe the 
real low-lying modes, and we meet an internal difficulty 
of theory. 

The author is grateful to Yu. V. Shtanov for the dis¬ 
cussion and the indication of an error in the analysis of 
the 5-function. I also thank the referees for helpful re¬ 
marks. 


IX. APPENDIX. COMPARISON OF THE 
SOLUTIONS FOR POINT AND ALMOST POINT 
POTENTIALS 


At a finite a, we seek a solution in the form 


'I'(a:) 


6i cos fca:-|-62 sin fca:, a: G [—L/2, —a[, 
die-^^+ d2e^^, |a:| < a, ( 74 ) 

oi cos fca;-|-02 sinfca:, xS]a,L/2]. 


Relation (1751) yields E = = \/{7 — k? (for |fc| < 

Vu). The boundary conditions and the sewing condi¬ 
tions for 'I'(x) and 'k'(x) yield the equations 

61 cos (fcL/2) — 62 sin (/cL/2) = 0, (75) 


Oi cos (fcL/2)-I- 02 sin (fcL/2) = 0, (76) 


61 cosfco — 62 sinfco = 516"^“-I-^26 ( 77 ) 

bik sin ka + 62 ^ cos ka = —Kdie”® -I- KC?2e“”“, ( 78 ) 


oi cosfco-I-02 sinfco = die -I-^2, ( 79 ) 


— aik sin ka +a 2 k cos ka =—Kdie + Kd 2 e^°‘. (80) 

They have two “branches” of solutions. For o —>■ 0 and 
k"^ <C U, we have 

I) di = d 2 = ail2, bi = oi, and 62 = —02 = 
01 /tan (A:L/2). The value of oi can be found from the 
normalization condition, and k satisfies the equation 

tan (fcL/ 2 ) = —fc/c. (81) 

II) di = —d2 = —ka2/2K —>■ 0, 5i = —oi = 0, 62 = 02 , 
02 is determined from the normalization condition, and 
k satisfies the equation 


The interatomic potentials U (r) are usually have a high 
repulsive barrier in the region r < 2 A and a shallow pit 
in the region 2 A < r < 5 A and tend asymptotically 
to zero, as r increases H- Is it possible to model a 
nonpoint high barrier with the 5-function? To answer, 
we compare the solutions for the wave functions and the 
energies obtained for both potentials. 

Here, we consider the simplest one-particle task: a par¬ 
ticle in the one-dimensional potential well —L /2 < x < 
L/2 (i.e., zero BCs at x = EL/2) with the potential bar¬ 
rier 

C/(|x|)= = (72) 

at the well center. Here, a <C L/2. The Fourier transform 
of such potential is J^(fc) = {2 RqUo sin ak)/ak. Bypassing 
to the limit a —>■ 0 , we have v{k) —>■ 2i?ot7o = 2 c = const, 
which corresponds to C/(x) = 2c5(x). Let us compare the 
solutions for almost point (arbitrarily small, but finite a) 
and point interactions. 

The Schrodinger equation reads 

- 4'"(x) -b C/(x)4'(x) = E^{x). (73) 


tan {kL/2) = 0, k = 2 t:1/L, Z = 1, 2,3,... (82) 

For the point potential U{x) = 2c5(x), we seek the 
solution of Eq. (1731) in the form (l74l) without the second 
row. We possess the BCs 

4/(±L/2) = 0 , 

the condition of continuity of 'L(x) at the point of the 
barrier 

4'(x = -5)|5 ^o = '^{x = 5)|5^o, 
and the equation 

[4''(x = 5) - 4''(x = -5)]= 2c4'(0), (83) 

obtained by the integration of the Schrodinger equation 
dZSl) on the interval x S [~5,5] (a similar equation arises 
also for N point bosons Q). These equations have two 
branches of the solutions: 1) E = k'^, bi = ai, 62 = 
—02 = oi/tan(fcL/2), and Eq. (ISTl) for k; 2) E = k"^, 
61 = oi = 0, 62 = 02 , and Eq. (1551) for k. They coincide 
with solutions (I) and (H) for almost point particles. We 
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note that though the function S{x) does not act on odd 
functions, such functions can be eigenfunctions of the 
Hamiltonian with the d-function. 

Let us consider the properties of solutions. For series 
(I), the lower level corresponds to the WF without nodes. 
The next levels correspond to the WFs with two, four, etc 
nodes. For series (II), the lower level corresponds to the 
WF with a single node. For the next levels, the WFs have 
three, five, etc nodes. By the theorem of nodes [3l|, the 
ground state corresponds to the WF without nodes, the 


first excited state to the WF with one node, the second 
excited state to the WF with two nodes, etc. Solutions 
(I) and (II) correspond to the theorem of nodes. 

Note that the eigenvalue E = oi the Schrddinger 
equation coincides with the value of ('I'o|IL|'I'o) both for 
the almost point interaction and for the point one. In 
the proof, it is necessary to consider that, for the point 
potential, ^'□(a:) has a discontinuity at the point x = 0 
(see (1551)1. 
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